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Effect of turbulent diffusion on iron abundance profiles 
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ABSTRACT 

We compare the observed peaked iron abundance profiles for a small sample of groups 
and clusters with the predictions of a simple model involving the metal ejection from 
the brightest galaxy and the subsequent diffusion of metals by stochastic gas motions. 
Extending the analysis of Rebusco et al. (2005) we found that for 5 out of 8 objects 
in the sample an effective diffusion coefficient of the order of 10 29 cm 2 s _1 is needed. 
For AWM4, Centaurus and AWM7 the results are different suggesting substantial 
intermittence in the process of metal spreading across the cluster. There is no obvious 
dependence of the diffusion coefficient on the mass of the system. 

We also estimated the characteristic velocities and the spatial scales of the gas mo- 
tions needed to balance the cooling losses by the dissipation of the same gas motions. 
A comparison of the derived spatial scales and the sizes of observed radio bubbles 
inflated in the ICM by a central active galactic nucleus (AGN) suggests that the 
AGN/ICM interaction makes an important (if not a dominant) contribution to the 
gas motions in the cluster cores. 
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1 INTRODUCTION 

X-ray spectroscopy is widely used to determine the metal- 
licity of the hot gas in galaxy clusters and groups. Using 
the high energy resolution of ASCA, BeppoSAX, Chandra 
and XMM-Newton the radial profiles of Fe, S, Si, Ca and 
other elements (e.g Buote et al. 2003b, Sanderson et al. 2003, 
Sun et al. 2003, Matsushita, Finoguenov & Bohringer 2002, 
Gastaldello& Molendi 2002, O'Sullivan et al. 2005, Fukazawa 
1994, Tamura et al. 2001) have been derived for a sample 
of objects. Outside the cluster core regions the metallicity 
of the intracluster medium (ICM) is on average one-third of 
the solar one and it does not seem to evolve up to a redshift 
of about 1 (e.g Mushotzky & Loewenstein 1997, Tozzi et al. 
2003). This suggests an early enrichment by Type II super- 
novae (SNII; e.g. Finoguenov et al. 2002). The metallicities 
in the core regions demonstrate much more diversity. For 
the so-called 'non-cooling flow' clusters the metallicity and 
the surface brightness do not vary strongly across the core 
region. For another group of so-called 'cooling flow' clusters, 
the metallicity and the surface brightness are both strongly 
peaked at the center. The clusters from the latter group al- 
ways have a very bright galaxy (BG) dwelling in their cen- 
ters, which makes these galaxies a prime candidate for pro- 
ducing the peaked abundance profiles. These peaked profiles 
were likely formed well after the cluster/group was assem- 



bled (e.g. Bohringer et al. 2004; De Grandi et al. 2004). The 
relative abundances of different elements indicate that Type 
la supernovae (SNIa) explosions and stellar winds within the 
BG have played a major role in the formation of the central 
abundance excesses (e.g. Renzini et al. 1993, Finoguenov 
et al. 2002). If the metals in the ICM are due to the galaxy 
stars then one would expect the abundance profiles to follow 
the BG light profile. However the observed metal profiles are 
much less steep than the light profiles, suggesting that there 
must be a mixing of the injected metals, a process which 
may help to diffuse them to larger radii (e.g. Churazov et 
al. 2003, Rebusco et al. 2005, Chandran 2005). Stochastic 
gas motions could provide such a mechanism of spreading 
metals through the ICM, provided that the characteristic 
velocities and the spatial scales are of the right order. Fur- 
thermore, the dissipation of the kinetic energy of the same 
gas motions could be an important source of energy for the 
rapidly cooling gas in the "cooling flow" clusters. 

In Rebusco et al. 2005 we estimated the parameters 
(velocity and length scale) of stochastic gas motions in the 
core of Perseus cluster (A426) required to spread the metals 
ejected from the central galaxy. We found that a diffusion 
coefficient of the order of 2 x 10 29 cm 2 s _1 is needed to ex- 
plain the observed abundance profile and that for turbulent 
velocities of about 300 km s _1 and eddies sizes of about 20 
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Name 


kT 


z 


NGC 5044 


1 


0.0082 


NGC 1550 


1.37 


0.0124 


M87 


2.2 


0.0044 


AWM4 


2.3 


0.0319 


Centaurus 


3.5 


0.0107 


AWM7 


3.7 


0.0176 


A1795 


6 


0.0639 


Perseus 


6.3 


0.0179 



Table 1. The sample: (1) name of the cluster , (2) mean ICM 
temperature in keV, (3) redshift. 



kpc the dissipation of turbulent motions compensates for the 
gas cooling losses. We now extend the same model to a small 
sample of clusters and groups. The sources we analyzed have 
mean temperatures which vary by a factor of ~ 6. 

The structure of the paper is the following. In section [5] 
we describe the sample of clusters and groups; in section |H] 
the structure and the ingredients of the stochastic diffusion 
model are described. The results are discussed in section [I] 
The last section summarizes our findings. 
We adopt a Hubble constant of Ho = 70 km s~ Mpc~ , 
M = 0.3 and ft A = 0.7. 



2 THE SAMPLE 

We have selected a small sample of nearby "cooling flow" 
clusters and groups (z < 0.07) having a sufficiently detailed 
information on the abundance distributions in the core re- 
gions. To test for any obvious trend with the mass of the 
cluster, we selected the objects having the mean tempera- 
ture ranging from 1 to 6 keV (see Table 0. 

2.1 NGC 5044 Group 

The NGC 5044 group of galaxies is relatively cool and loose. 
It is formed by a central luminous giant elliptical surrounded 
by about 160 galaxies, ~ 80% of which are dwarfs (Ferguson 
& Sandage 1990). There is evidence for multiple temperature 
gas components, that coexist within ~ 30 kpc. It is likely 
that the central heating is due to the presence of a central 
AGN (e.g. Buote et al. 2003). The X-ray emission is quite 
symmetric, indicating that there have not been recent strong 
mergers (David et al. 1994). 

2.2 NGC 1550 Group 

This luminous group is more relaxed than other bright low 
temperature ones (e.g. NGC 5044). The BG is a lenticular 
galaxy and it is almost isolated in the optical band (Garcia 
1993). 

Sun et al. (2003) suggest that NGC 1550 entropy profile 
shows signs of nongravitational heating (e.g. a recent out- 
burst). Moreover they point out the large role of the central 
galaxy in affecting the surrounding gas temperature. 
If this is the case, then the same nongravitational processes 
may help to distribute the metals within the group core. 



2.3 M87 Galaxy 

M87 is the dominant central galaxy in the Virgo Cluster. 
Because M87 is very close and bright it was used to study the 
role of SNIa and SNII in enriching the ICM (e.g. Gastaldello 
& Molendi 2002, Matsushita et al. 2003, Finoguenov et al. 
2002). 

There is clear dynamical evidence of the presence of a 
supermassive black hole in its core (e.g. Macchetto et al. 
1997), with a one-sided jet (e.g. Schreier et al.1982, Owen et 
al. 1989, Sparks et al. 1996). At least two gas components 
with different temperatures are observed: the hotter compo- 
nent is almost symmetric around the central galaxy (Ma 
tsushita et al. 2002), while the cooler component forms ex- 
tended structures correlated with the radio lobes. There are 
clear signs of an AGN/ICM interaction in this source both 
through the shocks and sound waves (Forman et al. 2005, 
2006) and through mechanical entrainment of the cool gas 
by the bubbles of relativistic plasma (Churazov et al. 2001). 

2.4 AWM4 Cluster 

This poor relaxed cluster consists of 28 galaxies centered at 
the dominant elliptical NGC 6051 (Koranyi & Geller 2002). 
The cluster's stellar component is comparable with that of 
M87 and NGC 5044, but the iron abundance profile is natter 
than in these sources and there is no evidence for multiphase 
gas (O'Sullivan et al. 2005). The density is not strongly 
peaked at the center and no strong temperature drop is 
seen. While it is not a prominent cooling flow cluster, it 
does contain a single cD galaxy and we include this cluster 
to see what will be the outcome of the analysis applied to 
this object. 

The central galaxy hosts a powerful AGN (4C +24.36) 
which could be driving the gas motions in this cluster. 
O'Sullivan et al. (2005) assert that in AWM4 there might be 
global motions of the galaxies in the plane of the sky: these 
galaxy motions could contribute in spreading the metals. 

2.5 Centaurus Cluster 

Centaurus is the third nearest bright cluster (after Perseus 
and Virgo). The asymmetric X-ray structure around the cen- 
tral cD galaxy (NGC 4696) gives evidence for dynamical 
activity (e.g. Allen & Fabian 1994). This cluster is very in- 
teresting because it houses the most prominent abundance 
peak known (e.g. Fabian et al. 2005, Sanders & Fabian 2002, 
Ikebe et al. 1999). The analysis of the cool Ha filaments in 
its core indicates previous episodes of radio activity (Craw- 
ford et al. 2005). Fabian et al. 2005 and Graham et al. 2006 
have recently calculated the effective diffusion coefficient in 
the picture of stochastic turbulence: we discuss our results 
in comparison to these earlier works in 14.21 

2.6 AWM7 Cluster 

AWM7 is a poor bright cluster, part of the Perseus-Pisces 
super-cluster and centered about the dominant elliptical 
galaxy NGC 1129 (e.g. Neumann & Bohringer 1995). It is 
elliptically elongated in the east-west direction and it follows 
the general orientation of the Perseus-Pisces chain of galax- 
ies (Neumann & Bohringer 1995). The 30 kpc shift of the 
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X-ray peak from the optical peak indicates that AWM7 has a 
cD cluster in early stage of evolution (Furusho et al. 2003). 
The determination of the effective radius of NGC 1129 is 
quite controversial (see Table 2). This uncertainty does not 
affect our estimates, however. 

2.7 A1795 Cluster 

This compact and rich cluster is believed to being dyna 
mically relaxed (e.g.Tamura et al.2001). On the other hand 
Ettori et al.(2002) suggest that the core has an unrelaxed 
nature, consistent with the detected motion of the central 
cD galaxy MCG +5-33-5 (Oegerle & Hill 1994, Fabian et 
al. 1994, Markevitch et al. 1998). There is a strong cooling 
flow in this cluster with a substantial temperature drop and 
a strong density peak in the core (e.g. Briel & Henry 1996, 
Allen et al. 2001). 

2.8 Perseus Cluster 

Perseus (A 426) is the brightest nearby X-ray cluster and it 
is one of the best-studied cases of cool core clusters, together 
with M87 and Centaurus. Its central elliptical galaxy (NGC 
1275) dominates in the optical light up to a distance of ~ 
100 kpc. It hosts a moderately powerful radio source, 3C 84 
(Pedlar et al. 1990). In the core region there is a complex 
substructure in the temperature and surface brightness di 
stributions including depressions in surface brightness due to 
rising bubbles of relativistic plasma (Churazov et al., 2000, 
Fabian et al. 2000) and quasi-spherical ripples (Fabian et 
al. 2003a, 2006). Optical Ha filaments, whose origin is not 
well understood, seem to be drawn up by the rising bubbles 
(Fabian et al. 2003b, Hatch et al. 2005). 



3 THE MODEL 

The model is essentially the same used in Rebusco et al. 
(2005). For completeness we reproduce its main features be- 
low. 

3.1 Diffusion of metals due to stochastic gas 
motions 

We assume a static gas distribution in the gravitational po- 
tential of the cluster. The gas density and temperature di 
stributions are known and assumed to be not evolving with 
time. We then suppose that the gas is involved in stochas- 
tic motions which do not affect the density and temperature 
distributions on time scales much longer than the time scales 
associated with the gas motions, but only spread the metals 
through the ICM. Due to such motions, the metals injected 
in the center of the cluster will be spread out off the center. 
We treat this process in the diffusion approximation: 

^ = V • (DnVa) + S, (1) 
at 

where n = n(r) is the gas density, a = a(r, i) is the iron 
abundance and S = S(r, t) is the source term due to the 
iron injection from the BG. D is the diffusion coefficient, of 
the order of ~ y , with v being the characteristic velocity of 
the stochastic gas motions and I their characteristic length 



scale. Once the diffusion coefficient is specified, eq^can be 
integrated. In what follows we considered constant diffusion 
coefficients: the effect of having a diffusion coefficient as a 
function of radius is discussed in Rebusco et al. (2005). 



3.2 Adopted light, gas density and metal 
abundance profiles 

We list the parameters of the adopted profiles in Table 2: we 
used the existing analytical approximations (when available) 
or made our own fits to the original data (all the references 
used are given in the Table). 

The light distribution of the central cD galaxies is mod- 
eled here by a simple Hernquist profile (Hernquist 1990): 
this is an acceptable approximation for the purpose of this 
study. For each source we compared the light distribution of 
the BG with the light distribution due to all the other gala 
xies excluding the BG. The brightest galaxy dominates up 
to a distance of ~ 100-150 kpc, hence in the calculations we 
assume (unless explicitly stated otherwise) that the central 
excess of metals is produced by the central galaxy alone. 

The iron abundance profile is approximated with a sim- 
ple /3-profile: a(r) = a(0) [l + (r/r a ) 2 \ , where a® is the 
solar abundance (Anders & Grevesse 1989). The adopted 
functional form should work only in the central region, where 
the abundance excess is present. This form also neglects 
completely the central abundance "hole" observed in some 
sources (e.g. Schmidt et al. 2002, Bohringer et al. 2001). The 
nature of this abundance hole remains unexplained and it 
may be due to a not adequate modeling of the emission from 
the very central region rather than due to a real decrease of 
the metal abundance. 

For the electron density profile we used a single (5 model 
again: n e (r) = n(0) [l + {r/r c ) 2 ] 3/ ' 2 ' 3 . The hydrogen num- 
ber density is assumed to be related to the electron number 
density as nn = n e /1.2. In some cases (marked with an 
asterisk) a double /3-profile is reported in the original pa- 
pers: here we use the parameters of the component which 
provides the dominant contribution for the range of radii of 
interest. Only in the case of Centaurus the double profile 
have been used, because both the components are impor- 
tant. For Perseus cluster we adopted the same profiles as in 
Rebusco et al. (2005). 



3.3 Iron enrichment 

Much of the metals in the cluster ICM were produced at 
early times by SNII explosions. They are likely to be evenly 
distributed through the bulk of the ICM producing a uni- 
formly enriched gas. The level of this uniform enrichment 
is uncertain - on average an iron abundance in the range 
0.2-0.4 of the the solar value is reported. In the subsequent 
calculations we subtracted this "abundance basis" in or- 
der to single out only the central abundance excess, which is 
believed to be formed later by the metal ejections from the 
BG alone (Matsushita et al. 2003, Bohringer et al. 2004, De 
Grandi et al. 2004). The assumed value of a,b is one of the 
main uncertainties in our model - for low values of the 
central abundance excess becomes more extended and the 
mass of iron in the central excess grows substantially. 
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Ram-pressure stripping can also contribute to the en- 
richment of the ICM over a long period of cluster evolution. 
According to recent simulations (Domainko et al. 2006) the 
cluster centers are enriched more strongly than the outer 
parts. It is not clear though how these results will be af- 
fected by the presence of cool and dense cores. Below we 
assume that ram pressure stripping produces broader abun- 
dance peaks than considered here and that the central abun- 
dance excess is solely due to the central galaxy. 

We assume that the central abundance excess is due to 
the type la supernovae and the stellar mass loss by stars in 
the BG. Both channels inject gas enriched with heavy ele- 
ments, in particular iron. Following Bohringer et al. (2004), 
the rates of iron injection by SNIa (eq. 2) and stellar mass 
loss (eq. 3) are: 



dM F 
dt 



( dM Fe 
\ dt 



Sli x 1.0 f |§ | , h . 



(2) 



= /,'(/) 0.105 > 1.0 ( ^§ ) \J. ,,i 



= 7Fe x 2.5 x 10~ n ( — 



-1.3 



(3) 



if) !/ "' 



where SR is the present SNIa rate in SNU (1 SNU - rate 
of supernova explosions, corresponding to one SN event per 
century per galaxy with a blue luminosity of 10 10 Lq), vfb = 
O.7M is the iron yield per SNIa, -y Fe = 2.8 x 10~ 3 is the 
mean iron mass fraction in the stellar winds of an evolved 
stellar population, tH is the Hubble time. The expression for 
the stellar mass loss was adopted from Ciotti et al. (1991), 
assuming a galactic age of 10 Gyr. The stellar mass loss 
contribution to the hydrogen content of the ICM can be 
neglected, as its effect is important only in the central ~ 
10 — 20 kpc (see Fig. 1). The time dependent factor R(t) = 
(t/tn)^ k ^ takes into account an increased SNIa rate in the 
past (Renzini et al. 1993), with the index k ranging from 
1.1 up to 2. We assume a fiducial value for the present day 
SNIa rate of 0.15 SNU (Cappellaro, Evans, & Turatto 1999). 
Hence the total iron injection rate within a given radius r 
can be written as 

fdMpe\ _ ( L B (< r) 
I dt 



SN 



T B 



(4) 



The procedure of evaluating the diffusion coefficient in 
equation (1) was as follows. First the value of at was spe 
cified and subtracted from the observed abundance profile. 
The iron mass in the remaining abundance excess was then 
calculated and compared with the total amount of iron pro- 
duced by the central galaxy during its evolution from some 
initial time to the present (by integrating equation 4 over 
time). This provides the constraints on the parameters of 
the enrichment model (SR, k, t ag e) (see Rebusco et al. 2005 
for an extended discussion). For each source we explored 
diffe 

rent models, with a b € [0.2, 0.4], k G [1.1, 2] SR G [0.07, 0.34] 
and tage G [5,10]. In particular, we fixed a b — 0.2,0.3 
and 0.4; then for each abundance basis we found a sam- 
ple of enrichment parameters: e.g. by setting k — 1.1, 2 and 
SR — 0.15, 0.3, and calculating the t age necessary to get the 
right amount of iron excess, then by fixing SR and t age and 



calculating k, finally by calculating SR from four combina- 
tions of k and t age . A representative set of these models is 
listed in Table 4. 

Equation (1) was then integrated over time from tH — 
tage. till now starting from a zero initial abundance. The re- 
sulting profiles were visually compared with the observed 
abundance peaks and the diffusion coefficient was adjusted 
so that the predicted and observed profiles would agree rea- 
sonably well (see Fig.l). In this figure for each object in the 
sample the dashed line shows the expected iron abundance 
due to the ejection of metals from the galaxy if no mixing 
is present. In most cases such profiles are far too steep than 
the observed abundance excesses (shown with a solid line). 
The dotted line in Fig.l shows the expected profiles calcu- 
lated for the same parameters of iron injection model, but 
with the additional effect of diffusion. The actual values of 
the diffusion coefficient used in each case are listed in Table 
3. Clearly allowing for metal mixing leads to a reasonable 
agreement between the observed and the predicted abun- 
dance profiles. 

3.4 Cooling and heating 

As shown above the observed and predicted abundance pro- 
files can agree reasonably well if the metals are allowed to 
diffuse through the ICM. Assuming that the diffusion coef- 
ficient derived above is due to stochastic gas motions, one 
can cast it in a form D ~ C\vl, where C\ is a dimensionless 
constant of the order of unity, v is the characteristic veloc- 
ity of the gas motions and I their characteristic length scale. 
Thus the above analysis provides an estimate of the product 
of v and I for each object. 

3.4-1 Dissipation of turbulent motions 

Here we further assume that the dissipation of turbulent 
motions is the dominant source of energy for the rapidly 
cooling gas in the cluster cores. As in Rebusco et al. 2005 
we suppose that the dissipation rate is simply equal to the 
gas cooling rate Tdiss ~ r coo ;. The gas heating rate due 
to the dissipation of the kinetic energy can be written as 
Tdiss ~ C2PV 3 /I, where p is the gas density and C2 is a di- 
mensionless constant, while the cooling rate F cooi = n 2 A(T) 
can be estimated from the observed gas temperature and 
density at any given radius. 

This requirement of the cooling and heating losses pro- 
vides a second constraint on the combination of v and I. Thus 
fixing the diffusion coefficient and the heating rate one can 
easily estimate v and I, that we intend as the characteristic 
sizes and velocities of quasi-continua intermittent vortexes. 

Below we use the expressions for D and Fdiss from Den- 
nis & Chandran (2005) (see references therein): 



r^is 



D = Civl ~ Q.llvl 

- C 2 pv 3 /l ~ 0.4pv 3 /l. 



(5) 
(6) 



To evaluate the heating rate for each object from the sample 
we use the gas parameters near ro = r coo i/2 (where r coo i is 
the radius at which the cooling time is comparable with 
the Hubble time). We then used D from Table 3 and set 
1 diss n A(T), evaluated at r . 
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Figure 1. Comparison of the observed and the expected iron abundance profiles for each source: the set of enrichment parameters used 
are listed in Table 4. The solid line shows the abundance profile adopted in this paper, from which a constant value of a b is subtracted. 
We assume that this central abundance excess is mainly due to the metal ejection of the central galaxy. For comparison we show the 
expected iron abundance (short dashed line) due to the ejection of metals from the galaxy. The expected profile was calculated assuming 
that the ejected metal distribution follows the optical light. In all the objects the expected abundance profile is much more peaked than 
the observed profile (due to the contribution of the central galaxy) suggesting that some mechanism is needed to spread the metals. The 
dotted line shows the profile derived with the same parameters of iron injection, but with the additional effect of diffusion. The long 
dashed line corresponds to the maximum abundance (it can be obtained from eq. 4), beyond which our approximation of neglecting the 
gas injection is not valid anymore. 
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Figure 2. Range of the characteristic velocities v and spatial scales I of the gas motions which provide the necessary diffusion and 
dissipation rates. Along the thick dot-long dashed line the diffusion coefficient D = C\vl is constant and equal to the values listed in 
Table 3, while the thin dot-long dashed lines show the effect of varying C\ = 0.11 by factor of 1/3 and 3. Along the the dotted line 
the dissipation rate is equal to the cooling rate at ro = r coo [/2; the dot-short dashed line is for ro = r coo [ and the long dashed for 
7*0 = r coo [/4 . At the intersection of these bands, two conditions are satisfied: i) the gas cooling is balanced by the dissipation and ii) the 
diffusion coefficient is of the right order. The vertical solid lines show the sound speed in the gas c s for r = r coo i/2 and half of it. In the 
plots where only one line is present, it corresponds to c s /2. 
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This approach of choosing ro = r coo i/2 implies that 
the cooling rates used will not be drastically different from 
source to source (as it would happen if for instance one uses 
a fixed value of ro for all the objects in the sample) . The re- 
sulting constraints are shown in Fig. 2. The thick dot-dashed 
line shows the combinations of v and I which give the same 
diffusion coefficient, while the thin dot-dashed lines show 
the effect of varying C\ by factor of 1/3 and 3. The effect 
of varying Ci by factor of 1/3 and 3 was estimated in Re- 
busco et al. (2005). Along the the dotted line the dissipation 
rate is equal to the cooling rate at ro = r coo ;/2; the dot- 
short dashed line is for ro = r coo i and the long dashed for 
ro = r coo ( / 4 . The intersection of the two bands gives the lo- 
cus of the combinations of I and v such that on one hand the 
diffusion coefficient is approximately equal to the required 
value (Table 3) and on the other hand the dissipation rate 
is approximately equal to the cooling rate at r coo i/2. 

3.4-2 Turbulent mixing vs turbulent dissipation 

Since cluster atmospheres are characterized by positive gra- 
dient of the specific entropy, stochastic motions should also 
lead to the heat flow into the central region (e.g. Dennis & 
Chandran, 2004). Following their work we estimate the rate 
of heating due to the turbulent transport of high-specific- 
cntropy gas into low-specific-entropy regions: 



r tt = V ■ (DpT gas X7e), (7) 

where T gas is the temperature and e = Cyln^/jp^ -1 ) is 
the specific entropy (adiabatic index 7 = 5/3, specific heat 
at constant volume CV = Sk/2 nn/p )• 

First we consider the case of a diffusion coefficient inde- 
pendent of the radius. A typical example of Ftt radial depen- 
dence (NGC 5044), for the diffusion coefficients taken from 
Table 3, is shown in Fig. [3] While the heat flow is always 
towards the center, the volume heating rate Ftt is positive 
inside a certain radius and it becomes negative (i.e. cooling) 
in the outer regions. This is of course an expected result - for 
constant diffusion coefficient the heat flux r 2 DpT gas ^ ~ rp 
(sec eq0 and there is likely to be a radius near which this 
function is almost independent from the radius. As it is seen 
in Fig. E]the turbulent transport (dotted line) calculated for 
D — 9 x 10 28 cm 2 s" 1 is an important source of heat in the 
innermost 10 kpc, but Ftt quickly declines with the radius. 
Obviously the diffusion coefficient has to increase with the 
radius if the turbulent heat transport F tt balances the cool- 
ing over a broad range of radii. For comparison the heating 
rate due to the dissipation of turbulent motions Tdiss (for 
the parameters taken from Table 3) is shown by the dashed 
line in Fig. [3] Note that the rate of dissipation is a very 
strong function of v and variations in v by a factor of 2 (for 
a fixed diffusion coefficient cx I x v) cause variations of Tdiss 
by a factor of 16 (see Fig. [3J. 

Chandran (2005) considered convective motions in the 
cluster cores driven by the cosmic rays injected by a cen- 
tral AGN. In his two-fluid (thermal-plasma and cosmic-ray) 
mixing-length theory approach, the characteristic length 
and velocity scales are functions of the radius. Assuming 
that the turbulent heat transport is the dominant source of 
heat we calculated the diffusion coefficient Do(r) (needed 
to balance heating and cooling at every radius) by integrat- 
ing the equation T t t = F coo i over the radius. The resulting 
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Figure 3. NGC 5044: the solid line shows the cooling rate as a 
function of radius; the dotted line is the heating rate Ftt evaluated 
for D = 9 X 10 28 cm 2 s _1 : the thin line corresponds to the part 
where Ftt is positive, while the thick one is for Ftt < 0. The re- 
maining lines represent the dissipation heating rate Fd i3s for the 
same D, but different velocity and length scales: I = l re f = 11 
kpc and v = v. re f = 245 km/s (short-dashed), I = l re f/l.8 and 
v = 1.8 v re f (dot-dashed), I = 1.8 l re f and v = v re t /1.8(dot- 
dashed). For the other sources, these plots show a similar quali- 
tative behavior. 



diffusion coefficient linearly increases with the radius and it 
reaches the value of D — 9 x 10 28 cm 2 s _1 at r ~ 5 kpc. 
Such radially dependent diffusion coefficient produces pro- 
gressively more and more efficient mixing with increasing 
distance and it leads to an abundance profile which drops 
to ~0.1 already at ~10 kpc (see Fig. 4). 

We next assumed that Fdiss ~ r coo ! at every radius 
and we calculated the quantity v :i /l as a function of the 
radius. We then considered several possibilities for the radial 
dependence of v and I and we derived the corresponding 
diffusion coefficients. 



Di(r) = civl[r), v — const 
Aa(r) = civ(r)l, I = const 
D 3 (r)= av(r)l(r), l(r)=0.3r. 



(8) 



The choice of l(r) oc r for Ds(r) is motivated by the 
work of Chandran (2005). For the first two models the con- 
stant v and / were fixed at the values taken from Table 3. 
The abundance profiles obtained by integrating equation 1 
using these diffusion coefficients are shown in Fig. [1] The 
case of Z?2(r) (i.e. constant length scale I and variable ve- 
locity scale v) fits the observed profile the best, while the 
other versions of D(r) yield profiles which strongly deviate 
from the observed one. 

From these plots one can conclude that it is possible 
to construct a model with radially dependent velocity and 
length scales such that the abundance profiles are approx- 
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Name 


r eff /kpc 


Ltot/Lq 


r a /kpc 


o(O)/o 


b 


r c /kpc 





n(0)10- 3 /cm 3 


NGC 5044 


10.05 [1] 


4.510 10 


[1] 


58.2 [2] 


1.07 


1.19 


15 [3] 


0.49 


10.66 


NGC 1550 


9.1 [4] 


3.510 10 


[5] 


0.8 [4] 


0.95 


0.14 


2.5 [4] 


0.35 


7 


M87* 


7.8 [5] 


5.310 10 


[21] 


20 [7] [8] 


0.85 


0.38 


2 [6] 


0.42 


130 


AWM4 


20 [9] 


5.510 10 


[10] 


94 


0.78 [11] 


0.89 


76 [3] 


0.62 


3.52 


Centaur us* 


13 [9] 


510 10 [22] 


25 [13] 


1.31 


0.48 


7.7[14] 


0.57 


165 


AWMT 


42 [9] 16.6 [12] 


1.410 11 


[15] 


7.2 [16] 


1.6 


0.3 


3.5 [17] 


0.25 


165 


A1795 


34[18] 


210 11 


2.7 [19] 


0.51 


0.21 


17 [19] [20] 


0.4 


50 


Perseus 


15 


1.610' 


.1 


48 


0.54 


0.18 


/ 


/ 


/ 



Table 2. Column (1) Name of the object (2)-(3) Effective radius and total blue luminosity of the central galaxy (4)-(5)-(6) Parameters 
of the iron abundance /3-profilc (7)-(8)-(9)Parameters of the gas density /3-profile (see the text for the precise form). For the sources 
marked with ( * ) the best fit was a double /3-model; here we report the parameters of the dominant one. References: [l]Buote et al. 
2004[2]Buote et al. 2003[3] Sanderson et al. 2003[4] Sun et al. 2003[5] dc Vaucoulcurs ct al. 1991 [6] Matsushita ct al. 2002[7] Matsushita 
et al. 2003[8] Gastaldello & Molendi 2002[9] Schombert 1987[10] Finoguenov ct al. 2001[11] O'Sullivan ct al. 2005[12] Bacon ct al. 1985 
[13] Fukazawa 1994[14] Ikebe et al. 1999[15] Peletier et al. 1990[16] Furusho et al. 2003[17] Neumann & Bohringer 1995[18] Schombert 
1988[19] Tamura ct al. 2001 [20] Ettori et al. 2002 [21] Arnaud ct al. 1992[22] Laine et al. 2003. 



imately reproduced and that the dissipation of turbulent 
motions compensates the gas cooling over a wide range of 
radii. In the simplest model of this type (D2 (r) above) the 
length scale is independent of the radius and the velocity 
scale is a slowly decreasing function of the radius. 



4 RESULTS AND THEIR DISCUSSION 

In Table 3 the estimated values of the diffusion coefficient 
D, the characteristic velocity v and the spatial length scale I 
are listed. The set of the parameters for the iron enrichment 
models (listed in Table 4, described in the Subsection 3.3) 
are chosen in order to get a reasonable approximation of the 
observed profiles (see Fig. 1): a factor of 3 smaller (or larger) 
diffusion coefficients produce too peaked (or too shallow) 
abundance profiles, inconsistent with the data. 

For each source we made a consistency check evaluating 
the characteristic value of the diffusion coefficient by com- 
paring the half mass radius of the observed central excess 
and the excesses produced by ejection and diffusion in the 
model. For small diffusion coefficients the metal distribution 
should follow the light distribution of the central galaxy (i.e. 
have the same effective radius), while for the increasing diffu- 
sion coefficient the effective radius should also increase. The 
values derived this way were consistent with those given in 
Table 3. 

We then used the observed abundance peaks in each ob- 
ject as initial conditions and we calculated the subsequent 
evolution of the abundance profiles by fixing D at the value 
from Table 3. In each case 3 times smaller (or larger) dif- 
fusion coefficient causes a quick steepening (or flattening) 
of the abundance peaks on time scales of a Gyr, while for 
diffusion coefficients of the order of those listed in Table 3 
the evolution is very slow on these time scales. 

The diffusion coefficients listed in Table 3 are mostly 
of the order of 10 29 cm 2 s _1 , apart from those obtained 
for AWM4, Centaurus and AWM7 (see next subsections). 
These values can be considered as an an upper limit on the 
effective diffusion coefficient set by the random gas motions, 
since other processes may contribute to spreading the metals 
through the ICM. 

Note that our simple model predicts the abundance to 



Name 


r coot /2 kpc 


D cm 2 s 1 


I kpc 


v km s 1 


V/C a 


NGC 5044 


20 


9 x 10 28 


11 


245 


0.51 


NGC 1550 


30 


1.2 x 10 29 


17 


219 


0.35 


M87 


20 


8 x 10 28 


10 


250 


0.34 


AWM4 


70 


4.5 x 10 29 


35 


378 


0.45 


Centaurus 


25 


2 x 10 28 


2 


250 


0.30 


AWM7 


50 


6 x 10 27 


1 


193 


0.18 


A1795 


50 


1.3 x 10 29 


14 


278 


0.25 


Perseus 


45 


2.4 x 10 29 


17 


410 


0.44 



Table 3. Estimate of the turbulence parameters for each source 
in our sample: name (1), half the cooling radius (2), diffusion co- 
efficient (3), length scale (4) and velocity (5-6) of the gas motion. 



be the highest in the center of the cluster. Thus the abun- 
dance "hole" observed in several well studied clusters (e.g. 
Bohringer et al. 2001, Schmidt et al. 2002) contradicts this 
conclusion. Potentially there are several explanations of the 
abundance hole. The simplest one supposes that the spec- 
tral models used for the abundance determination are not 
adequate for the central region, where a complicated mix- 
ture of different phases is present (e.g. Buote 2000). The 
resonant scattering and the diffusion of the line photons to 
larger radii is another possible explanation (e.g. Gilfanov et 
al. 1987, Molendi et al. 1998) . However the model advocated 
here assumes that the gas in the center is not at rest, but it 
is involved in complex pattern of motions, which should re- 
duce the effective optical depth (e.g. Churazov et al. 2004). 
Velocities of order the of few hundred km/s should make 
most of the lines optically thin. Less straightforward inter- 
pretations of the abundance hole are also possible, but they 
are beyond the scope of this paper. Future high energy res- 
olution observations with micro-calorimeters should be able 
to resolve this question. 

The sizes of the turbulent eddies obtained in this work 
are of the order of 10 kpc and the characteristic velocities 
of the order of few 100 km/s. In our model the velocity 
scale is determined from the combination of two quantities 
D ~ vl and T dtss ~ v :i /l: v ~ (D V di3S ) (1/i) . As a result 
the estimates of v come out very similar for all the objects, 
while I spans over a larger interval, compensating for all the 
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Figure 4. NGC 5044: in both plots the solid line represents the observed abundance profile from which the constant value a;, = O.3a0 
is subtracted. The dotted line is for D = Do(r), the short-dashed for D = Di(r), the long-dashed for D = D2(r) and the dottcd-dashed 
for D = D^(r)(sec subsection 3.4.2 for an explanation): they show the expected abundances obtained according to equation 1, starting 
from a fiat zero profile (left panel, T age = 5 Gyr) and starting from the present-day profile (right panel, after 4 Gyr from now). For the 
other sources, these plots display a similar qualitative behavior. 

peculiarities. We stress that these estimates have to be taken 
with caution, due to the simplicity of the model and the to 
uncertainties in the enrichment process. 



4.1 AWM4 

This object is not a cooling flow cluster, unlike the other ob- 
jects in the sample. It does contain however a central bright 
galaxy (NGC 6051) which should enrich the ICM with met- 
als in the same way as in cooling flow clusters. The enrich- 
ment will be even more efficient in AWM4 because the gas 
density is relatively small and even a modest amount of met- 
als should produce a prominent abundance peak. As one can 
see from Table 3 the diffusion coefficient required to explain 
the lack of such abundance peak in AWM4 is larger than for 
any other objects in the sample. The length scale is also the 
largest in the sample. This can be easily understood since in 
this source a large value of the diffusion coefficient is needed 
along with a modest heating rate. As it is clear from eq.(5) 
and (6) this can be most naturally done by setting a large 
value of the gas motions spatial scale I. 

If we would like to stay within the frames of our model 
then this is the object with the strongest level of gas mixing. 
There is a powerful radio source in NGC 6051 which clearly 
interacts with the ICM (e.g. O'Sullivan et al. 2005). It is 
not obvious however if the central AGN alone is responsible 
for the gas mixing and for the appearance of AWM4 as a 
non-cooling flow cluster. The alternative would be a recent 
merger followed by an intensive mixing of the ICM, but this 
case seems to be excluded by observations (O'Sullivan et al. 
2005). 



4.2 Centaurus 

The Centaurus cluster has a central abundance peak which 
is difficult to reproduce with the enrichment model presented 
above. In order to produce such a large amount of iron in 
the central excess of this cluster, a long accumulation time 
(of the order of 10 Gyr) and a high SNIa explosion rate are 
necessary (see Table 3). With these parameters we get a 
diffusion coefficient of ~ 1.5 x 10 28 cm 2 s _1 . 

A previous estimate by Fabian et al. (2005) obtains that 
the diffusion coefficient in Centaurus is lower than 6 x 10 28 
cm 2 s _1 , consistently with our findings. Our result is also 
consistent (within the uncertainties of the profiles) with the 
value found by Graham et al. (2005). Since Graham et al. 
(2005) use higher blue luminosity for the central galaxy, 
their choice of enrichment parameters is less extreme than 
our. They also adopted higher value of the central abun- 
dance peak (ao ~ 2<i0 ), derived from fitting the spectra 
with a cooling flow model, while our lower central abun- 
dance (ao ~ 1.3a©) is based on the two-temperature fits by 
Sanders & Fabian (2002). For their choice of parameters su- 
personic velocities are required to satisfy the equations (5) 
and (6). For lower luminosity (see Table 2), subsonic solu- 
tions are possible. However the results for Centaurus cluster 
remain quite controversial. 

One plausible explanation of the extremely steep abun- 
dance gradient in Centaurus is the recent disturbance of the 
gas in the cluster core. "Cold fronts" are observed in many 
clusters with cool cores (e.g. Markevitch et al. 2000) . These 
features are believed to form when the gas sloshing in the 
potential well brings layers of gas with different entropies in 
close contact: in a similar fashion layers of gas with differ- 
ent abundances can form a transient feature characterized 
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by a steep abundance gradient. In the frame of the model 
discussed here any steep abundance gradient translates into 
strong limits on the diffusion coefficient. 

4.3 AWM7 

Taking the results of our analysis at face value, this clus- 
ter does not require strong stochastic motions to match the 
observed and predicted abundance profiles. There is also no 
evidence for a strong radio source in it(e.g. Furusho et al. 
2003). From this point of view it is tempting to conclude 
that AWM7 is now passively cooling and that the activity 
of the central AGN is only about to start (in another cool- 
ing time ~7x I0 8 years). Then our estimate of the velocity 
and the length scale based on the assumption of a balance of 
the gas cooling would overestimate the combination v 3 /l and 
underestimate I. Note however that Furusho et al. (2003) re- 
ported an unusual substructure in the metal distribution in 
the AWM7 core, which is at odds with the assumption that 
gas motions are almost absent in this cluster. Given that in 
several clusters an abundance hole is observed it is possible 
that the substructure in AWM7 has a similar nature. 



4.4 Radio bubbles and stochastic gas motions 

Vazza et al. (2006) simulated non-cooling flow clusters, fin 
ding that the energy due to turbulent bulk motions scales 
with the mass of the cluster: similar calculations had been 
done by Bryan & Norman (1998). In these simulations the 
turbulence is a result of the cluster formation, while in the 
core of CF clusters one can expect that the AGN is playing 
a dominant role in driving the turbulence: as a consequence 
the mass scaling of the two processes may be different. 

Allen et al. (2006) found recently a tight correlation 
between the estimates of the Bondi accretion rate on to the 
central black hole and the jet power (estimated from the 
observed expanding bubbles) in X-ray luminous elliptical 
galaxies (similar calculations have been previously done by 
Di Matteo et al. 2001, Churazov et al. 2002 for M87 and 
by Taylor et al. 2006 for NGC 4696). At the same time 
the comparison of the jet powers and the cooling rates (e.g. 
Birzan et al. 2005) also shows a reasonable agreement. These 
results suggest that the amount of energy provided by jets 
(in the from of bubbles of relativistic plasma) is regulated 
by the gas properties in such a way that the gas cooling 
is approximately compensated. One can compare the rough 
estimates of the turbulence length scales (of the order of 
few - few tens kpc , see Table 3) with the observed sizes 
of radio bubbles in cluster cores (e.g. Dunn & Fabian 2004, 
Dunn, Fabian & Taylor 2005, Churazov et al. 2001). Broadly 
the numbers are of the same order. E.g. the observed radio 
bubbles sizes in M87 vary from ~ 1 to ~ 10 kpc (depending 
if one refers to the inner lobes or to the torus-like eastern 
or the outer southern cavities). In the Perseus, A1795 and 
Centaurus clusters the observed bubbles sizes also fall in the 
range few - few tens kpc. The characteristic velocities of the 
bubbles motions are of the order of ~ 0.5c s and they are also 
broadly consistent with the values listed in Table 3. From 
this point of view the observed bubbles seem to have all the 
necessary characteristics required to drive the gas motions 
needed to spread the metals. 



5 CONCLUSIONS 

We show that the iron abundance profiles observed in a small 
sample of nearby clusters and groups with cool core are 
consistent with the assumption that metals are spreading 
through the ICM with effective diffusion coefficients of the 
order of 10 29 cm 2 s -1 . This value does not show any obvious 
trend with the mass or the ICM temperature of the system. 

For characteristic velocities of the order of 300 km/s and 
length scales of the order of 10 kpc, the dissipation of the 
gas motions would approximately balance the gas cooling. 
When observations are available, our estimates of the turbu- 
lent length scales are consistent with the observed buoyant 
bubbles size, pointing at the central supermassive black hole 
as the likely origin of the core gas motions. 

Two objects in our sample require very different diffu- 
sion coefficients: for AWM7 very small level of gas mixing is 
needed, while for AWM4 the mixing has to be very intense. 
It is possible that these two objects hint to the intermit- 
tence of the gas mixing and that they represent two different 
episodes of the cluster evolution: AWM7 is only weakly dis- 
turbed now (but it will likely be disturbed soon by the onset 
of AGN activity), while AWM4 has been already mixed re- 
cently. O'Sullivan et al. 2005, using the analogy with M87 
(where an abundance gradient is present in spite of clear 
signs of a radio jets interaction with the ICM) suggested 
that it is unlikely that AGN activity is responsible for the 
absence of the abundance peak in AWM4. One cannot how- 
ever completely exclude the possibility that in the past there 
was a period of truly violent AGN activity in AWM4, which 
not only removed the central abundance peak, but also de- 
stroyed the whole cooling flow. 
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Effective diffusion coefficient: the models marked with the asterisk have an uncertainty due to the choice of the criterion to fix D of a 
factor %/ X 2, while for the other models the uncertainty is smaller. 
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